Resonance and limit cycle in a noise driven Lorenz model 
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The effect of an external noise on the Lorenz model is investigated near the onset of convection 
and near the Hopf bifurcation. We show the existence of a diverging time scale near the onset of 
convection and a resonance near the Hopf bifurcation. Our calculation provides an understanding 
of the noise induced stabilization of the limit cycle that had been observed numerically. 
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The noise induced phenomena have an intense inter- 
est in several disciplines ranging from physics, to chem- 
istry, and to biology Physical systems are usually 
not isolated from their environment, and the environ- 
mental influences on the system often appear as fluctu- 
ations. It is clear that the fluctuations are an integral 
part of the evolution of physical, chemical and biologi- 
cal systems and must be understood if we are to accu- 
rately and quantitatively describe the world around us, 
particularly at the small scale. In the past few years, it 
has become clear that fluctuations can actually be used 
constructively, by us or by nature, to produce organised 
behavior that is not possible in the absence of noise. Ex- 
amples where noise leads to organised behavior include 
stochastic resonance noise induced phase transi- 

tions 0, 0] , noise induced pattern formation j^, 0] , and 
noise induced transport [l|, 0| . The constructive role of 
noise is only possible in the non-linear non-equilibrium 
systems and is entirely the result of the intricate inter- 
play of noise and non-linearity away from equilibrium. 

The effect of noise on dynamical systems include noise 
induced hopping between multiple stable attractors[lo[ 
[HI and noise induced stabilization of the Lorenz attrac- 
tor [l6| near the threshold of its formation. We will 
focus on the second aspect here and as an example of 
the effect of noise on hydrodynamic instability. We con- 
sider the Lorenz model with an external noise source. 
Zippelius et. al.[l3| first studied the correlations in this 
model. We will supplement their work by focusing spe- 
cially on the situations where the control parameter is in 
the vicinity of an instability. This will lead to a noise 
induced resonance near the Hopf bifurcation point. We 
then focus on the effect that the noise can have on the 
limit cycle beyond the Hopf bifurcation point. Ordinar- 
ily the limit cycle is unstable, the numerical work of Gao 
et. al. [1] shows it can be stabilized by noise and here 
we explicitly show how such a stabilization is possible by 
drawing on the technique of Mclaughlin and Martin[l3|. 

We study the statistical properties of the following 



noise driven Lorenz model. 



i = a{y - x) + r)i{t) 
y = rx — y ~ xz + ri2{t) 
z = —bz + xy + 'q^{t) 



(1) 



where, a is the Prandtl number, 6 is a geometric fac- 
tor, and the Rayleigh number r is the control parame- 
ter in this model. The parameter cr = 10 and 6 = 8/3 
are held fixed at their standard values, as was originally 
done by Lorenz. rji{t) is an external white noise source 
with zero mean and < r]i{t)r]j{t ) >= Sij2D6{t — t ), 
z,j = 1,2,3. _D is the noise strength. D = describes 
the unforced Lorenz model. When the noise is absent 
the Lorenz system [3, [3] shows remarkable change in 
behavior depending on the control parameter r. In the 
conduction range, r < 1, the trivial steady state solu- 
tion is stable and loses stability to the other two describ- 
ing steady convection through a bifurcation at r = 1. 
Thus, for r > 1, there is a pair of stable fixed points, 
(±-\/6(r — 1), ±^b{r — 1), r — 1) and these in turn lose 
stability at r = rr = cr((T + 6 + 3)/(<t - 6 - 1) = 24.74 
through a Hopf-bifurcation. For r > tt, no stable steady 
state solution exists and the system has a strange attrac- 
tor. Depending on initial conditions, the solution may 
settle down on any of these three attractors. The trajec- 
tories are non-periodic and wander around in the vicinity 
of strange set of attracting points for r > and turbu- 
lence sets in. 

Now we study the effect of noise on the Lorenz model 
in the conduction, convection and turbulent regime. We 
calculate the time dependent correlation functions and 
compare with the behavior of the unforced model. 

In frequency space equation ([Ij reads 
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{—iuj + a)x{Lu) — ay{Lu) = rii{uj) (2) 
—rx{uj) + {—iuj + l)y{Lu) + ^ x{uji)z{uj - ui) = r]2{ijj) 



{—iLU + b)z{Lu) 



LJ2 



Linearizing around the steady state x ^ y = z ~ Q 
i.e. neglecting the non- linear terms in ([2]), we can cal- 
culate the correlation function by solving the linearized 
equation. The correlation of the x variable reads 



Cxxil^) = < x{uo)x{-Uj) > 

(a;2 + i) + ^2 



(3) 



{[uj^ + [r - l)<j]^ + {a + l)^uj^) 



The time dependent correlation function C^x (t) can be 
written as 



dLue-"^^Cxx{uj) 



(4) 



and is the sum of two exponential function with two 
inverse relaxation times which are independent of the 
strength D of the fluctuating force. Clearly, the actual 
relaxation time of the non-linear system ([T]) does depend 
on the strength of the fluctuating force. The correlation 
function of x gets damped out faster with increasing noise 
strength D. The correlation functions Cyy{t) and Czz{t) 
follow the same behavior in the regime r < 1. 

Near, r = 1, the correlation time goes as 



r oc 



r - 1 



(5) 



i.e. the relaxation time becomes infinitely big as r = 1 
is approached- a sign of critical slowing down. 

Now, for r > 1, Fourier analysis of the Lorenz model 
breaks down in its present form. This is very similar 
to what happens when one enters a symmetry breaking 
phase in critical phenomena. Accordingly we need to go 
to the shifted variable ui,U2, defined as iti — x — xq, 
U2 = y -yo and = z - zq. 

The Lorenz equation takes the following form 

{-iuj + a)ui{u;) - au2{LLj) = rii{uj) (6) 
-ui{uj) + {-iuj + l)u2{cj) + XqU3{lj) = ?72(w) 
-yoUi{uj) - xoU2{uj) + {-iuj + b)u3{LL>) = 773(0;) 

Now, the correlation function takes the following form 

< ui{llj)ui{—u!) >= (7) 
i-LU^ + bry + {b+ lyco^ + a'^{9 + g;^) + a%{r - 1) 
w2[u;2 _ + cr)]2 + [cj2(cr + 1 + 5) _ 2(7&(r - 1)]2 



The time dependent correlation function can be cal- 
culated from Cmut(t) = (iwe~*"*C„i„i (w). Their 
time dependence is determined by the three poles in the 
complex frequency plane. One is purely imaginary re- 
flecting the exponential decay in the correlation func- 
tion, the other two have flnite real parts of the opposite 
sign, reflecting spiral motion around one of the attractors 
(±a;o, ±yoj ^o)- This gives rise to the oscillatory behav- 
ior of Cuxui (t) superimposed upon the exponential decay 
of the correlations caused by the crossing between stable 
fixed points. In the absence of random force, the system 
is attracted in general to one of the stable fixed points 
(±Xo,±yo, zq), depending on its initial condition. If the 
random force is applied, the trajectories are no longer 
confined to one of the steady state points. The exponen- 
tial decay shows the motion from one fixed point to the 
other. The oscillatory motion slows down the decay of 
correlation function i.e. the memory effect of the initial 
state. 

At the Hopf-bifurcation point, i.e. r — the 
correlation function takes the following form 



< ux{uS)ux{-ijj) >= (8) 
i-uj'' + br)^ + {b+ Ifio^ + a^jb^ + c^2) + a^^^ _ i) 

[l^2_^2)]2[^2 + (^+ 1+5)2)] 

where, ujo — b{(7 + ^^r^rr^)- Now, the real time cor- 
relation function Curui{t) = dwe~^'^*'Cuiux{<^) exists 
as a principle value. Hence, the real time correlation 
function goes as 



(9) 



This is a consequence of the fluctuating force. It is sim- 
ilar in appearance to the resonance in a simple harmonic 
oscillator subjected to a sinusoidal force. 

Approach to this resonance is of the following 



< MiUi >= lim ie*"* 
c->o e 



(10) 



when, r 



CT(b+a-+3) _ 

In the absence of noise there is no periodic state above 
the Hopf bifurcation point. This is because the hmit 
cycle is unstable. In this case however, a periodic state was 
observed when the stochastic force was turned on. We 
will try to understand this on the basis of perturbation 
theory in terms of small noise strength e. We return 
to Eq.(8). Writing r = ro + Are, xa = xoo + e, where 

cAr 
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The zeroth order solution is found from 
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FIG. 1: Time series for the Lorenz equations with r=24.72, 
D= 0.001 as obtained by Gao et. al. 
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(14) 



M,o = i?eyl,e™* + G^j{t - t')f^j{t')dt" (15) 



where, Ai = A, A2 = ^{a - iw), A3 = ;^[1 - 
(i-iwo){cr-iwo) j Q_ found in frequency space from 
the particular integral of Eq. fM]) . At 0{e), we have 



where the operator L is given by 





(16) 



L = 



d_ 

dt 



a -a 
^ + 1 a;oo 



-a;oo -xoo + ^ 



(12) 



Consistent ordering in powers of e, leads to the expan- 
sion 



U3 



e^/^ I "20 





e^^^ I "22 



(13) 



and the driving term in the above equation has contri- 
butions which are; 

(i) a time independent term, 

(ii) a time periodic term with frequency 2wo, 

(iii) a term which is a product of a periodic term of fre- 
quency Wo and a stochastic term rj and 

(iv) a term which is a product of two stochastic terms. 
The solutions for un is a contribution of these driving 

terms with suitable coefficients. 

Having obtained un , we can now go to 0(e3/2) 

we 

have 
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(17) 



The driving terms on the right hand side need to be anal- 
ysed. The second term is time periodic with frequency 
wq and has an amplitude proportional to {Ar)A. In the 
first term, we need to analyse a typical term (unuao e.g.) 
to understand the different time dependent components 
of the drive. The term U30 is periodic with frequency 
wq. Two of the components of un as shown above are 
constant in time and periodic with frequency 2wo. The 
product M30w.11 then has terms with frequency wq and the 
amplitude of these terms are proportional to A^ . The op- 
erator L has a zero mode at frequency wq, consequently 
in order to allow a finite solution of Eq. (fT7| . the right 



hand side has to be orthogonal to the left eigenvector of 
L with zero eigenvalue. In the absence of the stochastic 
drive, this is the entire calculation and the result is that 
Ar — —0^A^, where /? is a constant. Thus there is no 
real value of A if Ar is positive and the limit cycle of the 
Lorenz model cannot be seen. 

In the presence of stochastic term, things change be- 
cause Mio (Ea. (|14p ) has a term proportional to rj and un 
has the two terms which are listed in (ii) and (iv) above. 
If we consider UioUji then there are two terms whose solu- 
tion is Arjrje™°* . Since we are in the presence of stochas- 
tic terms we can only talk about averages over 77. Conse- 



quently, we first need to average Eq. (jl7p over rj and only 
terms involving product of two Ty's will survive. After the 
averaging the term with structure Aryrye™"* will acquire 
the structure ^e™"* and will be a part of the dangerous 
term on the right hand side. The orthogonality condi- 
tion now leads to an equation Ar + N ~ — /J^A^, where 
N is the extra contribution coming from the stochastic 
term. In this particular case N is negative and hence 
for Ar < N, we can see the limit cycle stabilized and 
this is the mechanism which allows the time series to be 
periodic in a region of r, where there is no stable limit 
cycle. 

In closing, we have studied the statistical properties of 
noise driven Lorenz model. We have seen critical slow- 
ing down at stationary bifurcation point. At r > 1, 
the Fourier analysis breaks down due to lack of time 
translational invariance in present form of Lorenz model. 
An interesting resonance appears at the Hopf-bifurcation 
point, which shows the the noise induced stability. We 
also analysed the noise induced stability of limit cycle in 
the Lorenz model. 
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